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Abstract 

We study scaling properties of the localized eigenstates of the random dimer 
model in which pairs of local site energies are assigned at random in a one 
dimensional disordered tight-binding model. We use both the transfer matrix 
method and the direct diagonalization of the Hamiltonian in order to find how 
the localization length of a finite sample scales to the localization length of the 
infinite system. We derive the scaling law for the localization length and show it 
to be related to scaling behavior typical of uncorrelated Band Random Matrix, 
Anderson and Lloyd models. 
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1 Introduction 



It is well known that in one-dimensional (ID) disordered models even small amount 
of disorder leads to an exponential localization of all eigenstates 0] . On the other 
hand, recent studies of quasi-lD polymers have shown that short-range correlations 
embedded in a random sequence can lead to appearance of fully transparent states 
0, ||]. In ref. j| in particular, various organic disordered systems were quoted with 
electrical properties. The prototypical case is that of the Random Dimer Model (RDM) 
0, |[] where (in the context of a tight-binding Hamiltonian) pairs of adjacent energy 
levels are assigned at random, leading to two-site correlations in an otherwise random 
model. 

Since for infinite samples fully delocalized states appear only for specific energy 
values, there is no Anderson transition in the usual sense (see also [H). However, the 
number of transparent states for finite samples was found to be proportional to square 
root of the length of the sample ||, |[]. This fact is related to the divergence of the 
localization length in infinite samples when the energy approaches some critical values 
H |7], ||. Therefore, these states may be important for conducting properties of finite 
samples || [T(| . 



In infinite samples the Anderson transition can be characterized in terms of the 
localization length; the latter is commonly defined from the decay of amplitude of 
eigenstates in the limit |n| — > oo, where n is the site label in the tight-binding picture. 
Contrary to what happens in infinite samples, the global properties of eigenstates of 
finite samples cannot be characterized in the same way; one needs to use other quanti- 
ties (such as the inverse participation ratio), that are valid both for finite and infinite 
samples. Then, through the use of scaling conjectures, one can link the properties of 
eigenstates in infinite samples to those in finite samples. In the theory of disordered 
solids, the scaling approach proved to be extremely useful in describing the conduc- 
tance and its fluctuations (see e.g. 0). A similar approach has been recently 
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used in the theory of quantum chaos when describing localized eigenstates random on 
a finite scale [IB], 113. Such eigenstates also arise in the quasi-lD models with random 



potential. Extensive numerical and analytical studies (see e.g. [ITR |TB| and references 



therein) have revealed remarkable scaling properties of eigenstates, which seem to be 
of quite generic nature. 

In this paper we study the RDM of finite size from the point of view of scaling 
properties of its eigenstates. The question of the relevance of the above mentioned 
results to models with correlated disorder is far from being trivial since short range 
correlations may cause significant difference in the structure of eigenstates, when com- 
pared with those for random potentials. In the next Section 2 we briefly describe the 
RDM and discuss different definitions of localization length, which are used in our 
numerical simulations. In Section 3 we present numerical data on scaling properties 
of eigenstates in the center of energy bands. In this case, the localization length in 
infinite samples has been obtained by the transfer matrix method. In Section 4, we 
study the energy region near the critical values E cr by making use both numerical and 
analytical treatment of the localization length. Finally, in Section 5 we give a short 
summary of our investigation. 

2 Finite size scaling approach to random dimer model 

Our starting point is the ID Schrodinger equation in the tight-binding approximation, 

z ' ~~jT~ = e nC n (t) + c n+ i{t) + c n _i(t) , (1) 

where c n (t) is the probability amplitude for an electron to be at site n and e n is the 
local site energy. By making the transformation c n (t) = exp(-iEt) x n one can obtain 
the equation 

E(f n = ip n+ l + e n ip n + ip n -i , (2) 

for the eigenvalue E and the corresponding eigenstate <p n (E) . In what follows, we 
consider the RDM |3j, |T^] which implies short range correlations in the sequence of 



random e n . In this model there are only two values of e n viz. or e# that appear 
in pairs in the sequence e n 's . In other words, in order to create the dimer chain the 
pairs AA and BB each with energies €a, e_B respectively are distributed at random. 
We take for simplicity the probabilities, of occurrence of the paired to be equal, i.e. 
Vaa = Vbb = 1/2 . 

The RDM has been well studied for an infinite chain (see, e.g. @, ||, [| [7], ||, |i7f). 



Its basic feature is that for two values of the energy E cr = €a or e# its eigenstates are 
extended. In the vicinity of these values and for ea,b = Ya — e s| which is less than 
the critical value eJ B = 2 , the localization length defined through the exponential 
decay of the amplitude of eigenfunction diverges as loo(E) ~ 1/E 2 (for the specific 
value €a,b = 2 , the singularity law looi^E) ~ 1/E holds instead, see details in |], |7|. |J). 
In spite of the fact that for other values of E inside the spectrum the localization 
length is finite, the influence of nearly-transparent states on the electronic properties 
of finite samples is strong (see e.g || [TtJ). This is related to the fact that the number of 
eigenstates with localization length larger than the size N of the sample is proportional 

to Vn. 

Our main interest is in the structure of eigenstates for finite samples, both in the 
center of energy band and in the vicinity of the critical energies E cr where localization 
length in infinite sample diverges. Unlike the more simple case of infinite samples, the 
meaning of localization length for finite samples is not clear. Below, we follow the 
approach developed in the theory of quasi- ID disordered solids which is based on the 



evaluation of multifractal localization lengths (see, e.g. [15]). One of the commonly 
used quantities in this approach is the so-called entropic localization length, defined 
through to the information entropy 7i of eigenstates, 



N 



H = -^2w n ln w n ; w n = \ifl\ (3) 



n=l 



where (p n is the n— th component of an eigenstate in a given finite basis. For eigenstates 
normalized as Xm w n — 1 ■> t ne simplest case of (p n = N' 1 / 2 results to the entropy 
equal to the maximum value, viz. 7i = ln(iV) . We define therefore the localization 



length l N as the number of basis states occupied by the eigenstate ip n ; the latter 
is equal to exp(7-^v) • Similar definitions have been used for the first time in ref. 
[ P"5|] where different characteristics of eigenstates have been discussed in solid state 
applications. One can see that in the other limit case of an exponentially localized 
state with cp n = Z^ 1//2 exp(— |n — no | /Zoo) , the quantity In is proportional to Zoo , viz. 
l N « eZoo (assuming Z^ -C N ). One should note that, generically the amplitudes (p n 
fluctuates strongly with n and thus the coefficient of proportionality between In and 
Zoo depends on the type of fluctuations. 

To study the properties of chaotic states, localized on some scale in the finite 
basis in [|13| , [L4j] it was suggested to normalize the localization length In in such a way 
that in the extreme case of fully extended states the quantity In is equal to the size of 
the basis N . In such an approach, the entropic localization length zjj is defined as 

1$ = Nexp(<H N > -H ref ) (4) 

In Eq. the ensemble average < .... > is introduced over the number of eigenstates 
of the same structure. The normalization factor TC re f has the meaning of an average 
entropy of the completely extended random eigenstates in the finite basis, therefore, it 
can be easily found analytically |l~3| , 



Ke/ = ^(y + l)-^(^ln(^) (5) 

where ip is the digamma function and the distribution of components <p n is assumed to 
correspond to the Gaussian Orthogonal Ensemble (GOE). From @ one can see that 
for N ^> 1 the entropic localization length of random eigenstates, defined simply as 
exp{H re f), is approximately 2.07 times less than iV ; this result is due to gaussian 
fluctuations in the components ip n . 

Analogously, the whole set of localization lengths 1$ can be defined in the fol- 
lowing way [15], |Tj|: 



l^=N(^^)^; q>2 (6) 

Kef 
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where 

N 

n = EW (7) 

n=l 

and PM is the average value of P q for the reference ensemble of completely extended 
states. One should note that for the particular case q = 2 the quantity P2 is known 
as the participation ratio; it is widely used in solid state physics. In the limit case 

(2) 

of the GOE, one can find that P re j. — 3/N , therefore, the inverse participation ratio 
(P 2 )~ 1 , which is commonly taken as the definition of localization length, for random 
eigenstates is 3 times as less as the "actual" length N . 

In fact, the above expressions for the localization lengths 1$ is defined through 
the 2q— th moments of a distribution of components (p n of eigenstates; non- normalized 
to pM quantities (^) are well known in the multifractal analysis of wave functions. 
Such normalization turns out to be extremely important when establishing scaling 
properties of eigenfunctions. Indeed, by normalizing the localization lengths 1$ to the 
size N of the sample, 



P. 



(?) 



_ l _n_ 

q ~ N 



one can expect, in the spirit of renormalization theory, that the set of dimensionless 
parameters f3 q is the proper quantity to characterize generic properties of eigenstates for 
finite samples. According to the scaling conjecture in the modern theory of disordered 



solids, it was assumed |20[ that for quasi- ID disordered models described by Band 
Random Matrices the quantity (3 q depends on the scaling parameter A only, which 
is the ratio of the localization length for the infinite sample, to the size N of the 
sample itself. Therefore, the scaling relation can be written as 

A = /.(A); A = ^ (9) 

Detailed studies, both numerical and analytical, have confirmed this conjecture for 
different models like the Kicked Rotator Model and Band Random Matrices (see, e.g. 
14, [15, 1(| and references therein). Moreover, the scaling function f q (X) has also been 



found. 



Our main question is whether the relation of the type (^) is also valid for our 
dimer model with the correlated disorder. The first nontrivial question arises about 
the reference ensemble for the computation of the average entropy 7i re / . Indeed, in 
application to ID Anderson type models (see details in PT| ) the reference ensemble 
can not be chosen as an ensemble of full random matrices, like the GOE. This point is 
related to the fact that in the Anderson case fully extended states are not gaussian ran- 
dom functions but just plane waves which arise for zero disorder. In the dimer model, 
the situation is even more complicated due to strong dependence of the localization 
length on the energy. However, and this is our expectation, in spite of the presence 
of the extended states at the critical energies, scaling properties of the eigenstates in 
the dimer model of finite size N are of generic type discovered for ID and quasi- ID 
disordered models. 



For this reason and in the spirit of Ref. [21, [EJ, we define the normalization 
factors H r ef and P r gl from the solution of Eq. (|2|) for the zero disorder, e n = , 



E k = 2cos J ^ J , (10) 



h / / 2 \ nkn , . 

^ = i\—i) am —v 

with k, n = 1,...,N. The entropy Ti. re f of the above eigenf unctions in the large N 
limit has the same value for every eigenvalue E k , i.e., 

H ref = ln(2A) - 1 , (12) 

and correspondingly, 

- 4 < 13 > 

3 Scaling properties of localization lengths in the 
center of energy bands 

Since all results depend on the difference €a,b — \ca ~ € b\ but not on the actual values 
tA and ts separately, we can set = for simplicity. One should stress that both 
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localization lengths l^ and iff are functions of the energy E . For this reason, in our 
numerical experiments we consider ensembles of states specified by the values of the 
energy E in a small window AE and by different realizations of random on-site energies 
e n . We choose the size of the energy window in such a way that for every chosen value 
of €b the localization length is approximately constant inside this window (in all 
the cases ^ < 0.06 ). 

too 

In order to study scaling properties of the localized eigenstates we have used the 
transfer matrix method for infinite chains as well as the direct diagonalization of the 
Hamiltonians that one associated with Eq. ([!]) , for finite chains of size N. To find the 
localization length l^ we have studied the asymptotic behavior of the random matrix 
product n M n , where M n is defined through the relation 

£ n+ i = M n U, M n =(^ V ™ ~^y, v n = E — e n (14) 

for the vector £ n = (x n , x n _i) with the matrix M n known as the transfer matrix. Then 
the localization length is the inverse Lyapunov exponent 7 ; the latter is evaluated 
as the exponential decay rate of an initial vector £ x , 

N 

1 In n \M n U 

1 -=i = £5Ln %\ ■ (15) 

Although the Lyapunov exponent 7 for finite N depends on a particular realization of 
the disorder, for N — > 00 it converges to its mean value. For the above calculations we 
have used samples of length 5 x 10 5 for relatively large values of ee and up to 4 x 10 6 
for small values of e# . 

To reveal scaling properties of localization length for finite samples, we have 
computed two localization lengths 1$ and iff according to the relations discussed in 

,(2) 



the previous Section, with the normalization factors TC re f and -P re j- in the forms (|T2| ) 
and (|T3|). In the computations of these lengths, the energy window was taken in the 
center of the spectrum, around the value E = for e# equals 2, 1.8, 1.6, 1.2, 

1, 0.8, 0.6, 0.4, 0.35 and for the fixed value of N. The width of the windows has 



been numerically chosen to provide a small change of localization length inside any of 
windows. The values of (3\ and (3% are obtained by the averaging over an ensemble of 
random samples of size N = 100 -i- 800 for the values of cited previously As a result, 
the total numbers of eigenstates in the energy windows were more than 1000. 

All the data have been fitted to the scaling function (3 q found for quasi-lD dis- 



ordered models 15 



In fact, this scaling relation is exact only for q = 2 , however, for other cases of small 
values q , including q = 1 , it is very close to the correct one (see details in [|15|] ). 

Numerical data reported in Fig.l give clear evidence of a scaling of the type fllB]) . 
The fitting parameters c q are equal to C\ = 2.80 and C2 = 1.55 . From this figure one 
can see that the behaviour of f3 q is very different in the two limits of very localized 
{(3 q <C 1) and extended {(3 q ~ 1) eigenstates. The dependence fllBj ) has the remarkable 
property which can be seen in other variables, 

H-M^r), x- H ±) (it) 

which are more convenient when considering the whole region of both very localized 
and extended states. Indeed, in these variables the scaling has extremely simple form, 

Y q = a q + b q X (18) 

with b q — 1 and a q = ln(c g ) . The data for the scaling in variables Y,X are presented 
in Fig.2. The fitting parameters b 12 are found to be quite close to 1 i.e. b\ = 1.02 and 
62 = 0.98, for this reason in the Fig.2 we put b\ = 62 = 1- The remarkable result is that 
the above simple scaling relation holds in a very large region of the scaling parameter 
A = loo/N . According to the fit to the dependence flIB|) , the values ai j2 are: ai = 1.05 
and ct2 = 0.45, which gives Aa 12 = a± — a<i = 0.6. It is very interesting that these 
values of a\£ are the same as for common Anderson model [2~2"| in the center of the 
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energy band. This fact is very important in establishing the link between the RDM 
and Anderson models of finite size. 

It is of special interest to relate the entropy localization length 1$ and the local- 

(2) 

ization length l N associated with the inverse participation ratio. Their interdependence 
is shown in Fig. 3. We see that they are approximately equal for very localized and very 
extended states. It is also clear that fa is always less than fa since < P^ +1 \ due 
to the definition of Eq. fl7|). Using the definition of Eq. (|T6| ) one can find the relation 

between fa and fa: 

1 + (c - I) fa Ci 

4 Scaling of localization lengths near the critical 
energy. 



In the previous section we have shown that the scaling law (|ISD , found for fully 
disordered ID and quasi-lD models, also holds in our dimer model of finite size when 
considering eigenstates in the center of energy bands. In a sense, this property may 
be expected since far from critical energies where the localization length diverges, 
the eigenstates are assumed to be similar to that known for disordered models. The 
important question is whether this scaling holds for all energies inside the band, in 
particular, near the critical energies E cr = €a,£b- Direct numerical computation of 
the localization length through the transfer matrix method is very difficult in this 
energy region due to very weak convergence of Lyapunov exponents. For this reason, we 
have used the analytical expression which was derived for near the critical energies 
in an approach developed in || : 

Here, the factor Q stands for the probability for the pair e n = e n+ \ = eg to appear, 
and S is defined by the relation E = ts — 8 ~ e B |2B|. We remind that in our case 
Q = 1/2 and €a = has been assumed for the simplicity. From the above expression 
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one can find that if the value of e# is far from the stability border Eb = 2 , and 
the distance A = 2 — e# is large compared to 5 = es — E , the localization length 
diverges as 



2A 

QS 

In the other limit case of = 2 we have 



loo « 7T7^ 5 < A < 1 (21) 



loc ~ 7^, 5 < 1, A = (22) 



It is interesting to note that the same expressions (|2T| ) and (^) are obtained in Ref. 
[ O by assuming that localization length is determined by the reflection coefficient 
from a single pair e n = e n+ i , embedded in a perfect lead. It is of interest to check how 
accurate are estimates found in |§ and |17| if to apply them for any energy inside the 
band. 

To find the localization lengths and 1$ for finite samples of the size N, 
we have used the same approach described in the previous section, by examining the 
eigenstates with energies in a small energy window AE < lO~ 2 E cr near the critical 
energy E cr = e#. Yet, since in the region of critical energies loo{E) and thus the 
localization properties of eigenstates, depent from the energy in a singular way, (see 
Eq. (pip and Eq. (p2])) we took from the energy window only the eigenvector with the 
corresponding eigenstate which is closer to E cr (but always differant from it, E ^ E cr ). 
This was a natural choice in order to study statistical properties of eigenstates with the 
similar localization properties (i.e. the eigenstates just near the totally extended one). 
The average values of V$ have been obtained by statistical averaging over an ensemble 
of more than 3000 samples of the size N = 100 -j- 800 with different pair-correlated 
disorder. The results are reported in Fig. 4 together with the fit to the expression fllTf ). 



One can see a quite good scaling of the form ([16]), in spite of fluctuations which are 
much larger in this energy region compared to that in the center of bands. The fitting 
coefficients c\ = 2.20 and C2 = 1.06 are slightly less than those in the band center. 
This fact may be explained by an approximate character of the analytical expression 
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(one should also note that for the values of (3 q very close to the limit f3 q = 1 the 
computational errors are very large). 

In Fig. 5 the same data are represented in the variables (0), with the fit corre- 
spondent to the dependence (|i~8|) . It is interesting that in spite of a slight difference for 
the coefficients a\ = 0.75 and 02 = 0.08 in comparison for those found in the center of 
bands, the shift Aai^ = a\ — 0,2 = 0.67 remains almost the same (compared to 0.6). 

5 Summary 

We have studied a ID tight-binding model with binary on-site disorder that is randomly 
assigned in every two sites. For such a model we know that there exist two special 
energies E cr at which transparent states appear [Q, [7), §|. For other energies, but 
close to critical ones, the localization length is very large, leading to nearly-transparent 
states that are of great importance in the conducting properties of finite samples. This 
property is quite different from genuine disordered models of Anderson type. 

Our numerical study of random dimer models of finite size deals with the scaling 
properties of the eigenstates. This study was motivated by the remarkable scaling law 
that has been found for different ID and quasi-lD models, both dynamical (Kicked 



Rotator on a torus |13|, |Tj]) and disordered (ID Anderson and Lloyd models pl| , 



and quasi-lD models |15, 20]). These latter results indicate that eigenstates in finite 
samples with disorder have generic properties, regardless of the details of the disorder. 

The main result of our computations is that scaling properties of eigenstates of 
finite dimers are of the same type as for the disordered models mentioned above in spite 
of the existence of nearly-transparent states. In particular, both entropy localization 
length and localization length from inverse participation ratio normalized in the proper 



way, follow the universal scaling law of Eq. ( jig) 



The scaling relation of Eq. (|16|) can be also represented in a very intriguing form 
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Ill 

+ T7^-^7 (23) 



l%\e,E) M\e,E) 1^(0, E) 
which still has no physical ground. In Eq. (^), 1$ (e,E) represents the localization 
length for a finite sample and finite disorder, l^(e,E) is the localization length for 
an infinite sample with the same disorder and iffi®, E) is the localization length for a 
finite sample with zero disorder. One should stress that all three localization lengths 
are defined in the same way given through expressions of Eq. (^) and (0) . One can see 
that the form Eq. (|23"D is parameter independent; the same form holds also for the ID 



Anderson and Lloyd models (see [21], |2"2fl). 

In our numerical study the energy window has been chosen in the middle of 
the spectrum as well as close to the critical energy, giving the same scaling form 
(|16|). The slight difference in the coefficients c q for these two energy regions seems 
to indicate that the analytical expression ([H]) needs some correction related to an 
additional dependence on the energy when the latter is not close enough to the critical 
one. Our results indicate that the same scaling is expected to hold for other values 
of the energy inside the band. One should note that the scaling flip] ) (or, the same, 
(|23"D) can be used to check the accuracy of expressions for the localization length 
in dependence of the parameters E and , if for some values of these parameters the 
scaling function [3 q is found with a high accuracy. 

It is of quite interest to check the scaling behaviour of localization lengths corre- 
sponding to the higher moments q > 3 in @. Analytical treatment [|15[ for disordered 
models have shown that the scaling law (|1(J) approximately holds also for higher mo- 
ments. The correct expression for (3 q (X) is known only in the limit case of very localized 
(A < 1) and extended (A > 1) states. It has the same form ([18]) with b q = 1 but with 
different values of a q in these limits (see details in [filj). On the base of our results 
for q — 1, 2, it is quite natural to expect that for the dimer model the correspondence 
to the analytical predictions |TjJ also hold for higher moments, however, this question 



remains open. 
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Finally, we should comment that the results obtained in the present work can 
generalized to cases with correlation blocks larger than dimers, viz. m-blocks with 
m = 3,4,5, ... In these more general cases the following surprizing result holds: Given 
an arbitrary distribution of correlated blocks with even length, i.e. an arbitrary distri- 
bution of dimers, quatromers, sextomers, octamers, etc, with the same energy e, that 
populate a lattice with sites that have some other energy value, there is always a reso- 
nant energy E cr = e that corresponds to a delocalized state. This result can be easily 
deducted from the general expressions of ref. |§. On physical grounds, we expect the 
localization properties of the eigenstates of this system to follow similar scaling laws 
to the ones derived in the present work. 
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Figure Captions 



Fig.l Scaling of f3 q as a function of the localization ratio A = loo/N for the RDM 
with Q = 0.5 and = 0. The energies E are taken in an energy window centered 
at E = e B /2 for values of e B = 2.0; 1.8; 1.6; 1.2; 1.0; 0.8; 0.6; 0.4; 0.35. Smooth curves 
correspond to the dependence ( |IED with c q as a fitting parameter. 

(a) Scaling for /3i with c\ = 2.80 ; 

(b) Scaling for /3 2 with c 2 = 1.55. 

Fig. 2 The scaling of fti, /3 2 as a function of A = loo/N in the variables Yi )2 and X (see 
([17])) for the same values of the parameters as in Fig. 1. Straight lines (1) and (2) corre- 
spond to the expression (|i~8|) with a\ = 1.05; b\ — 1 and a 2 = 0.45; 6 2 = 1 respectively. 

Fig. 3 A plot of /3 2 as a function of It is interesting to note that the fitting curve 
has the same form ( p!9|) as the ones for the case of /3x i2 plotted in Fig.l as a function of 
A = loo/N. 

Fig. 4 The same as in Fig. 1 for the energies close to the critical one (E cr = eb) , for 
q = 1 (Fig. 4a) and q = 2 (Fig.4b). The values of e B are taken as e B = 1.8; 1.6; 1.4; 1.2; 
1.0. 

Fig. 5 The same as in Fig. 2 for the parameters of Fig.4. Straight lines (1) and (2) corre- 
spond to the expression (|T8|) with a\ = 0.75; b\ — 1 and a 2 = 0.08; b 2 = 1 respectively. 
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